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Atomic ensembles have shown to be a promising candidate for implementations of quantum 
information processing by many recently-discovered schemes. All these schemes are based on the 
interaction between optical beams and atomic ensembles. For description of these interactions, one 
assumed either a cavity-QED model or a one-dimensional light propagation model, which is still 
04 ' inadequate for a full prediction and understanding of most of the current experimental efforts which 

are actually taken in the three-dimensional free space. Here, we propose a perturbative theory 
to describe the three-dimensional effects in interaction between atomic ensembles and free-space 
light with a level configuration important for several applications. The calculations reveal some 
significant effects which are not known before from the other approaches, such as the inherent 
mode-mismatching noise and the optimal mode-matching conditions. The three-dimensional theory 
confirms the collective enhancement of the signal-to-noise ratio which is believed to be one of the 
main advantage of the ensemble-based quantum information processing schemes, however, it also 
shows that this enhancement need to be understood in a more subtle way with an appropriate mode 
matching method. 
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Recently, many interesting schemes have been proposed which use atomic ensembles with a large number of identical 
atoms as the basic system for quantum state engineering and for quantum information processing. For instance, one 
can use atomic ensembles for generation of substantial spin squeezing and continuous variable entanglement 

H-[|, for storage of quantum light j7|-[lT|, for realization of scalable long-distance quantum communication Jl2| , 
and for efficient preparation of many-party entanglement [ fj"3"| . The experimental candidates of atomic ensembles 
can be either some cold or ultracold atoms in a trap PJlC|], or a cloud of room-temperature atomic gas contained 
in a glass cell with coated walls d|],[Il|]. The schemes based on atomic ensembles have some special advantages 
compared with the quantum information schemes based on the control of single particles: firstly, laser manipulation 
of atomic ensembles without separate addressing of individual atoms is normally easier than the coherent control of 
single particles; Secondly, and more importantly, atomic ensembles with suitable level configurations could have some 
kinds of collectively enhanced coupling to certain optical mode (called the signal light mode) due to the many-atom 
interference effects. Thanks to this enhanced coupling, we can obtain collective enhancement of the signal-to-noise ratio 
(that is, the ratio between the controllable coherent interaction and the uncontrollable noisy interactions) compared 
with the single-particle case if we choose to manipulate the appropriate atomic and optical modes. This collective 
enhancement plays an important role in all the recent schemes based on the atomic ensembles. 

To describe the interaction between atomic ensembles and optical beams, in particular, to understand the collective 
enhancement of the signal-to-noise ratio, normally one assumed a simple cavity-QED model [[7PJl^,p| or a one- 
dimensional light propagation model ||,^| with an independent spontaneous emission rate for each atom. In contrast 
to this, most of the current experiments are done with free-space atomic ensembles coupling directly to the three- 
dimensional optical beams jj^J^JnJ . It is not obvious that the predictions from the simple models will always be valid 
for this real much more complicate experimental situations. To fully predict and understand the real experiments, 
one needs to answer various questions associated with the three-dimensional interaction effects, for instance, what is 
the inherent mode structure of the signal light when we have a definite geometry of the atomic ensemble? Can we 
achieve a good mode matching with the matching efficiency in principle arbitrarily close to one? What is the noise 
magnitude associated with density fluctuation (induced by the random initial distribution of the atom positions and 
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the random atomic motion) of the atomic ensemble? In the three-dimensional configuration, is there still the collective 
enhancement of the signal-to-noise ratio? 

In this paper, we will try to provide answers to the above important questions for a Raman type A-level configuration 
which is useful for scalable long-distance quantum communication |p"2"| and for many-party entanglement generation 
p3[ . In general, it is very challenging to build a full quantum theory from the first principle to describe the interaction 
between the many-atom ensemble and the infinite-mode optical field and to give definite answers to the above-listed 
questions. In this paper, we will use a perturbative approach by assuming that the Raman pumping laser is very weak. 
There are several motivations to use a perturbative approach: firstly, the schemes in jl2],[l3) for scalable quantum 
communication and for many-party entanglement generation can work in the weak-pumping limit, so the perturbative 
approach describes a useful realistic situation. Secondly, the perturbative calculations allow us to investigate the 
three-dimensional effects and to give definite answers to the questions listed above from the first principle without 
any doubtful approximation. Finally, we expect the three-dimensional effects should be at least to some extent 
independent of the power of the pumping laser, so the perturbative calculations could give at least some indications 
for the three-dimensional interaction picture in other regions. 

The calculations reveal some significant results which are unexpected from the simple models, for instance, it turns 
out that due to the density fluctuation of the atomic ensemble, there will be two sources of noise for the light-atomic- 
ensemble interaction: one is the spontaneous emission loss and the other is the inherent mode mismatching noise 
(here, by "inherent", we mean that this mode mismatching is not from any technical imperfection). Thus, we need 
to use two quantities to describe the signal-to-noise ratios, and need to keep a balance between these two sources of 
noise by appropriate mode matching methods to optimize the setup for applications. The intuitive mode matching 
method will result in a quite large inherent mode mismatching noise. It is better to use some other mode matching 
methods which reduce the mode matching noise at the price of increasing the spontaneous emission loss. This can 
optimize the setup for some applications, such as the ones in Refs. |l2] , |l3| ], since there the spontaneous emission 
loss has a far less important influence. The calculations in this paper also demonstrate that in the realistic three- 
dimensional configuration, one can still obtain large collective enhancement of the signal-to-noise ratio for atomic 
ensembles compared with the single-atom case, which is an important feature of this kind of systems. 

It is also helpful to make a comparison between the collective enhancement of the signal-to-noise ratio we consider 
here and some other collective optical effects, such as the superradiance. The similarity lies in that both phenomena 
involve many-atom interference. However, there are also important differences. For superradiance, it is a stimulated 
emission effect, and the enhancement is in the emission speed. At the initial stage, coherence is built between 
different atoms from the spontaneous emissions. After the coherence has been built, superradiance can be understood 
even from a classical interference picture. For the light-atomic-ensemble interaction considered here, we are still 
completely in the spontaneous emission region which needs a full quantum description, and the enhancement is not 
in the emission speed, but in the signal-to-noise ratio when we only manipulate and measure a definite atomic mode 
(called the symmetric collective atomic mode) and a definite signal light mode (which is the optical mode collinear 
with the Raman pumping light). The coupling between the above atomic and the optical modes are coherent for 
different atoms (which means for this coupling there is a certain phase relation between different atoms), while the 
coupling of these two modes to other atomic and other optical modes (called the noise modes) are inherent with 
random phase relations between different atoms (the randomness comes from the random atom positions). We have 
collective enhancement of the signal-to-noise ratio since the coherent coupling interferes constructively for different 
atoms. Of course, in the three-dimensional free-space, there are infinite optical modes, and the signal light mode 
changes continuously to other noisy optical modes when we continuously vary the solid angle. Due to the continuous 
change from the coherent coupling to the inherent coupling, we have non-trivial mode matching problem, and we get 
some inherent mode mismatching noise. 

We also would like to mention that there are some early important works on the transverse effects of the Raman 
interaction by using the three-dimensional light propagation equations [pr[ . However, it is not clear to us how to use 
this approach to describe the collective enhancement of the signal-to-noise ratio, and in particular, how to give definite 
answers to the questions listed above. There is also some unpublished effort to try to figure out the three-dimensional 
effects in another light-atom interaction configuration p|. The interaction picture is not yet clear to us. 

This paper is arranged as follows: in Sec. II, we explain the interaction scheme and the basic ideas of the applications 
of this interaction scheme for quantum information processing. With the applications in mind, we can focus our efforts 
on the most relevant quantities that we need to calculate. We will also explain in this section the basic interaction 
picture between the light and the atomic ensemble from our calculations. Then, in the next section, we will describe the 
theoretical model for this light-atom interaction from the first principle, and solve this model by using a perturbative 
approach. In Sec. IV, we will discuss the properties of this solution, in particular, we will define and calculate the 
mode structure of the signal light, the spontaneous emission inefficiency, and the inherent mode mismatching noise. 
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Sec. V is devoted to the discussion of the appropriate mode matching methods for real experiments. We show that 
by appropriate mode matching methods, we can keep a balance between the two sources of noise mentioned above in 
order to optimize this setup for applications. Finally, in Sec. VI we summarize the results. 



II. THE INTERACTION SCHEME AND ITS APPLICATIONS 



The basic element of our system is an atomic ensemble, which consists of a cloud of N a identical atoms with the 
relevant level structure shown in Fig. 1. A pair of metastable lower states \g) and \s) can correspond, for instance, to 
hyperfine or Zeeman sublevels of electronic ground state of alkali atoms. The relevant coherence between the levels \g) 
and |s) can be maintained for a long time, as has been demonstrated both experimentally and theoretically 

All the atoms are initially prepared in the ground state \g) through optical pumping. 




g 

FIG. 1. The relevant atomic level structure with a pair of ground or metastable states, and |e) , the excited state. 



The ensemble is then illuminated by a weak pumping laser pulse which couples the transition \g) — > |e) with a large 
detuning A, and we look at the spontaneous emission light from the transition |e) — > |s), whose polarization and/or 
frequency are assumed to be different to the pumping laser. There are two interaction configurations: in the first 
configuration, the pumping laser is shined on all the atoms so that each atom has an equal small probability to be 
excited into the state |s) through the Raman transition. In the second configuration, as is reported in the experiment 
||, the pumping laser is focused with a transverse area smaller than the transverse area of the atomic ensemble, so 
only part of the atoms arc illuminated by the laser at each instant. However, during the light-atom interaction period, 
all the atoms in the ensemble are moving fast, and they frequently enter and leave the interaction region. As a result, 
each atom still has an equal small probability to be excited into the state |s). 

As we will show in the next section, after the atomic gas interacts with a weak pumping laser, there will be a special 
atomic mode s s called the symmetric collective atomic mode, and a special optical spontaneous emission mode a s 
(which couples to the transition |e) — * |s)) called the signal light mode. The symmetric collective atomic mode s s and 
the signal light mode a s are defined respectively by 



No. 



(1) 



J /k°krf 3 k, 



(2) 



where au represents the the plane wave mode with the wave vector k (we have used the plane wave modes as the 
eigenmodes for the expansion of the optical spontaneous emission field). The operators <Zk satisfy the standard 
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commutation relations 



ak, a k , = S (k — k'), and is the normalized signal mode function whose explicit form will 

be specified in Sec. IV. We just need to mention here that a s represents the spontaneous emission light which is 
basically collinear with the pumping laser, with a distribution only over a very small solid angle. The particularity 
of the modes s s and a s comes from the fact that they are dominantly correlated with each other, which means, if 
a atom is excited to the symmetric collective mode s s , the accompanying spontaneous emission photon will most 
probably go to the signal light mode a s , and vice versa. There are still many other atomic modes in the ensemble and 
infinite other optical modes in the spontaneous emission field, and these atomic and optical modes are correlated with 
each other in a complicated way. However, the correlation between the good modes s s and a s is quite "pure", and 
these two modes are only weakly correlated with the other atomic and optical modes which contribute to noise. The 
correlation between the atomic mode s s and other optical modes contributes to the spontaneous emission loss, and 
the correlation between the signal mode a s and other atomic modes contributes to the inherent mode mismatching 
noise. The interaction picture for this system is schematically shown by Fig. 2. 




Other 
atomic modes 




Other 
optical modes 



FIG. 2. For a many-atom ensemble, The intuitive interaction picture between the atomic modes and the free-space optical 
modes, where cl represents the "good" correlation between the collective atomic mode and the signal light mode, c2 represents 
the spontaneous emission loss, and c3 represents the inherent mode mismatching noise. Broader connections stand for stronger 
correlations between the excitations in the corresponding modes. 



The applications of this system for quantum information processing exactly comes from the almost pure correlation 
between the modes s s and a s . If we neglect the weak correlations between the good modes s s , a s and the other noisy 
atomic and optical modes, that is, we neglect the two sources of noise illustrated in Fig. 2, all the noisy modes can be 
traced over, and we get effectively a two-mode problem. The excitations in the modes s s and a s can be both separately 
measured. For the signal mode a s , this is done through a single- photon detector with an appropriate mode- matching; 
and for the atomic mode s s , this can be done by first mapping the atomic excitation to an excitation in the signal 
mode through a repumping laser pulse p2| , and then detecting it again through a single-photon detector. Using the 
almost pure correlation between the modes s s and a s , one can generate some preliminarily entanglement between 
two distant atomic ensembles by only linear optics means, which forms the important first step for applications of 
this setup for different kinds of quantum information processing tasks detailed in Refs. Here, we will briefly 

explain the basic ideas on how to entangle atomic ensembles using the correlation between the modes s s and a s . This 
explanation helps us to define the most relevant quantities for the applications that we need to calculate. 

The setup for entanglement generation between the two distant atomic ensembles L and R is shown by Fig. 3. 
We apply simultaneously two short Raman pumping pulses on the ensembles L and R, respectively, so that for each 
ensemble the light scattered to the signal mode a s has a mean photon number much smaller than 1. The signal modes 
are then coupled to optical channels (such as fibers) through mode matching after the filters, which are polarization 
and frequency selective to filter the pumping light. The signal pulses after the transmission channels interfere at a 
50%-50% beam splitter BS, with the outputs detected respectively by two single-photon detectors Dl and D2. If 
either Dl or D2 registers one photon, the process is finished. Otherwise, we first apply a repumping laser pulse to the 
transition |s) — > |e) on the ensembles L and R to set the atoms back to the ground state \g), then the same Raman 
driving pulses as the first round are applied to the transition \g) — > |e) and we detect again the photon number of the 
signal modes after the beam splitter. This process is repeated until finally we have a click in either Dl or D2 detector. 
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FIG. 3. The schematic setup for generating entanglement between two distant atomic ensembles L and R. 



We show now how the entanglement is generated between the ensembles L and R if we successfully register one 
photon in Dl or D2 detector. To understand this, let us first look at one ensemble. Before the Raman pumping pulse, 
the collective atomic mode s s and the signal light mode a s are in vacuum states, which are denoted respectively by 
|0 a )j 1 p ). If we neglect the small noisy correlations, after the weak pumping, the state of the the modes s s and a s 
has the following form (see the next section for the detailed derivation) 



where p c -C 1 is the probability for one atom excited to the collective atomic mode s s , and o (p c ) represents the 
high-order terms with more excitations whose probabilities are equal to or smaller than p\. 

In Fig. 3, the pumping pulses excite both ensembles simultaneously, the whole system is thus described by the state 
\4>) L ®\4>) R , where L and \4>) R are given by Eq. (1) with all the operators and states distinguished by the subscript 
L or R, respectively. The two signal light modes are superposed at the beam splitter, and a photodetector click in 
either Dl or D2 measures respectively the operators a+a + or ala_ with a± = (a S L ± e t(p a sR ) /V2. Here, ip denotes 
the difference of the phase shifts in the two-side optical channels. Conditional on the detector click, we should apply 
a projection operator a+ or a_ onto the whole state \<p) L £g> \<p) R . If we neglect the high-order corrections in Eq. (3), 
the projected state of the ensembles L and R thus has the form 



which is maximally entangled in the excitation number basis. The generation of this kind of state forms the basis of 
further applications in quantum information processing J^,[D|. If we take into account the high-order terms in Eq. 
(3), the fidelity between the actually generated state and the ideal state (4) will decrease by an amount proportional 
to p c , and this decrease will contribute to the final fidelity imperfection of the schemes in Refs. |l^Jl^]. For the 
applications described by these schemes, we need to fix the fidelity imperfection to be small, which means that we 
should keep a small excitation probability p c by controlling the intensity of the Raman pumping laser. 

Now we look at the influence of the two types of noisy correlations illustrated in Fig. 2. The spontaneous emission 
loss can be quantified by the probability p sp0 n of the event that the scattered photon goes to some other optical modes 
instead of the signal mode a s while the accompanying atom is excited to the collective atomic mode s s . Without the 
spontaneous emission loss, for each round of Raman pumping, we succeed with a probability 2p c to get a detector 
click, which prepares the entangled state (4). However, the spontaneous emission loss means that even if an atom is 
excited to the mode s s (which has a probability p c for each ensemble), the accompanying photon, with a probability 
1 — Pspon, cannot be registered. For entanglement generation, the fidelity imperfection to the ideal state (4) is given 
by the excitation probability of the atomic mode s s , which is now p c j (1 — p S pon) instead of p c (we always use p c to 
denote the possibility of the good event described by Eq. (3)). To fix the fidelity imperfection to be small, in the 
presence of the spontaneous emission loss, we need to further decrease the excitation probability p c by a factor of 
1 — Pspon, which means that we have a smaller probability to register one signal photon. So, as a result of this noise, 
the preparation efficiency (the success probability of the scheme) is decreased by a factor of 1 — p sp on when we fix the 
fidelity imperfection. 

The inherent mode-mismatching noise can be quantified by the probability p mo dc of the event that an atom is 
excited to some other atomic modes instead of the collective mode s s while the accompanying photon goes to the 
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(4) 
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right signal mode a s . In this case, we can register a photon from the single-photon detectors, but the atomic mode s s 
will be actually still in the initial vacuum state |0 o ). So, this noise will add some vacuum component to the generated 
state between the ensembles L and R. In the presence of this noise, the generated state between the ensembles L and 
R is mixed with the form 



PLR (CO,?) = (c \0 a a ) LR (0 a a \ + \* V )? a 



(5) 



where the vacuum coefficient cq is basically given by the conditional probability p mo dc for this noise contribution. 
Actually, as has been shown in , there are other sources of noise which can contribute to the vacuum coefficient 
Co, such as the detector dark counts, and the detector inefficiency in the succeeding entanglement connection scheme 
detailed in Ref. fl^| . It has also been shown there the vacuum component noise will be finally automatically purified, 
and thus has no influence on the final communication fidelity. However, it has a more significant influence on the 
efficiency than the spontaneous emission loss. For applications in Refs. to get a better overall efficiency, it is 

better to keep a balance between the mode mismatching noise p mo de and the spontaneous emission inefficiency Pspon; 
with Pmode significantly smaller than p S pon- 

It is helpful to make a comparison here with the single-atom case. For a single atom interacting with the free-space 
light with the same level configuration as shown in Fig. 1, there is only one atomic mode given by s = \g) (s|, but 
there are still infinite optical modes. We can still identify the forward scattered optical mode a s as the signal mode. 
The interaction picture is then shown intuitively by Fig. 4. For a single atom, there is no inherent mode-mismatching 
noise, but the spontaneous emission inefficiency becomes much larger. If the atom is excited to the level |s) (that 
is, to the mode s) through the Raman laser pumping, the accompanying photon has a probability to go to all the 
possible directions jl6| , and thus has only a very small possibility to go to the signal mode a B . We will calculate in the 
following sections the spontaneous emission inefficiencies for the signal atom case as well as for the atomic ensemble 
case, and compare them with the different mode-matching methods to see whether there is collective enhancement of 
the signal-to-noise ratio for the many-atom ensemble. 



cl 
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FIG. 4. The intuitive interaction picture for the single-atom case with the corresponding symbols having the same meaning 
as Fig. 2. There is only one atomic mode, but the spontaneous emission loss represented by c2 becomes much larger compared 
with the good correlation cl. 



III. THE THEORETICAL MODEL FOR THE LIGHT-ATOMIC-ENSEMBLE INTERACTION AND ITS 

SOLUTION 

We will now go to the detailed description of the light-atomic-ensemble interaction. From this description, we 
will calculate the signal-light mode structure /k, the spontaneous emission inefficiency Pspon? 

and the inherent mode- 
mismatching noise p mo dc- In the interaction configuration shown by Fig. 1, the Raman pumping laser can be described 
classically by neglecting its small quantum variance, and is assumed to be propagating basically along the z direction. 
So, we can write its amplitude as e ge (r, t) — u (r, t) e l ( k oz-^ot) ^ wnere Wo = fc QC = 2ttc/Xo is the carrier frequency and 
u (r, t) is a slow varying function of the coordinate r and the time t. The field e se coupling to the transition |e) — > \s) 
should be described quantum mechanically, and we expand it into plane wave modes e se (r,f) oc J a k e i( - k - r - Ukt U 3 k 
with = |k| c, and at, the corresponding annihilation operator. Then, the Hamiltonian describing the light-atom 
interaction has the following form in the interaction picture (setting h — 1) 

H (t) = A £ <4 + g eg < g u (rut) e* k ° z (6) 



G 



N a 

^crl I g*aie- t[k - r '~ {uJ *-^ +u °° )t] d 3 k + H.c. 



i=l 



where the detuning A = Lo eg — luq with u eg = uj e — ui g , the frequency difference between the atomic levels, = (v\ 
(/i, v = g, e, s) are the transition operators of the «th atom, i\ is the coordinate of the ith atom, and g eg ,g^ e are the 
coupling coefficients which are proportional to the dipole moments of the corresponding transitions. The coefficient 
g^ e depends in general on the direction of the wave vector k by the dipole pattern. In Eq. (6), the carrier frequency 
of the spontaneous emission field e se (r, t) is given by loq — uj sg , and its relevant frequency width can be estimated by 
the natural width T of the excited level |e), which means that the modes ak with the frequency difference between 
Wk and too — u sg much larger than the width T will have negligible influence on the system dynamics. We neglect in 
Eq. (6) the spontaneous emission back to the level \g) since it is not important for our purpose (with the detection 
method specified in Sec. 11), and has no influence on all of our results. 

If both the natural width T and the frequency spreading of the pumping laser u (r, t) are significantly smaller than 
the detuning A, we can adiabatically eliminate the upper level |e) through the standard technique, and the resulting 
adiabatic Hamiltonian is given by 



|2 N a 



+H.c.}-^5>»>(r l ,t)| 2 , (7) 

i=l 

where Aa;k = ^k — ( w o — ^sg) an d Ak = k— /cqz with z , the unit vector in the z direction. We have neglected in 
Eq. (7) the Stark shift of the level |s) since it is much smaller than the other terms (The spontaneous emission field 
is much weaker than the pumping field). The last term in Eq. (7) (the Stark shift of the level |<?))can be eliminated 

if we make a phase rotation of the basis {|g), \s)} which will transform a\ g to a l sg e"^ 9ea ^ $o u ( r " T 'l dT _ Thus, we 

can simply drop off the last term, and replace u (r.j, t) by u' (r i; t) — u (r i: t) e^ 9 "^ /o l u ( r,,r -'l rfT . i n the following, 
for simplicity of the symbol, we will denote u' (r,, t) still by u (rj, t), and denote —g eg g^ e /A by a single coefficient #k- 
At the beginning, all the atoms and in the ground state \g), and all the optical modes ak are in the vacuum state. 
We denote this initial state of all the atomic and optical modes by |vac). Then, if the Raman pumping laser u (r, t) 
is short and weak enough, we can expand the final state |\&/) into a perturbative expansions. To the second order of 
the perturbation, the final state has the form 

to 



|* / ) = |l-ij^ H(r)dT 



It 

2 



f f 1 1 



|vac) (8) 



where io is the duration of the Raman pulse, and T [■ ■ •] denotes the time-ordered product. 

At this point, we should note that for an atomic vapor, the atom positions are randomly distributed, and during 
the light-atom interaction, the atoms are moving fast in the room-temperature p7| . Therefore, we should treat in 
the Hamiltonian (7) as stochastic variables. We make the following two assumptions about the properties of these 
variables: (i) For different atoms i and j, Vi and do not correlate with each other; (ii) Different stochastic variables 
Yi obey the same probability distribution which is determined by the geometry of the atomic cell. Note that these 
two assumptions are well satisfied by the room-temperature atomic gas. If the atom positions r.^ behave as classical 
random variables, we should take an average of the final state |^/) over the joint probability distribution of these 
variables, and the resulting state becomes mixed with 

Pf = (1%) <*/l> {r4} . (9) 

where the symbol (• • '){ r } denotes the average over all the variables r^. 

To take the average in Eq. (9), first we note that the Hamiltonian (7) can be written as 

H{t)=Y^H i {v u t) (10) 
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with 

Hi (rut) = o\ g J tfkajy^* { u (n, t) e- lAkr ' } d 3 k. (11) 

Only the term in the large bracket of the Hamiltonian Hi depends on the variable r^. By the property (i) of the 
variables r i7 we have 

(HiHj) {ri} = (Hi) {r . } (Hj) {r , } 

+ 6 v{( H ?){r i} -( H if{r i} }- (12) 

By the property (ii) of the variables r^, we know that (u (r^, t) e~ jAk ' ri )^- r y becomes independent of the atom index 
i, so the average of the Hamiltonian H (t) has the simple form 

(H(t)) {Ti} = VK Ss (13) 

x J 3k4e lA " kt < W (r 4 ,i)e- 4Ak ^) {ri} d 3 k, 

where s s is the symmetrical collective atomic operators defined in Eq. (1). We will write the initial vacuum state 
| vac) as a tensor product of the atomic part |vac) a and the optical part |vac) p , and denote the atomic state cr* s |vac) a 
simply by |s)j. With this notation, a combination of Eqs. (8)-(13) yields the following form for the averaged state pf 

Pf = (l-P2- Pc) |*eff> (*eff| + VlPn + O (p 2 ) , (14) 

where the effective pure state I'J'cff) (not-normalized) is 



*eff) = 1 - * / (H (r)) {ri} dr \vac), (15) 



with (vteffl ^eff) = 1 +Pc, and the noise component p 2 Pn is given by 



P2Pn = l^ 4 ( S \® J J dT l dT 2 J J d3k ' 

X3k.9 k -e t ( AwkT1 - A ^' T2 )aj c |vac) p (vac| a k , (16) 



x {(ufc.nK (ri^)e- < ( k - k ')' Pi ) {p } 

- (u (r u ri ) e - 4Ak ' r {ri} (« ('i. ^) e iAk '" ri ) {ri} } 



with the value of £»2 determined by the normalization of p n , i.e., by tr (p n ) = 1. In Eq. (14), 0(^2) represents the 
higher-order terms compared with the magnitude of p 2 and p c , which will be neglected in the following. We also 
ignored the second-order cross-terms proportional to algO^vac) (vac| or |vac) (vac| <J l gs a^ in writing Eq. (14), since 
they have no contributions to the quantities we will calculate in the below. The solution (14) serves as the starting 
point for discussions of various properties of this system in the following section. 



IV. PROPERTIES OF THE LIGHT-ATOMIC-ENSEMBLE INTERACTION 

Now, let us analyze the properties of the light-atomic-ensemble interaction revealed by the solution (14). First, we 
look at the effective pure state |* ff) by neglecting the noise component P2p n - The state \*S> e fi) has exactly the same 
form as the ideal state (3) which is the starting point for all the applications. To see this clearly, we can write \^ e s) 
as I'J'eff) = (l + \/Pc s t a t) l vac ); by defining the signal mode 

/ k 4d 3 k (17) 
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where the value of p c is determined by the normalization of the operator a\, i.e., by the condition [a s ,aj] = 1. 

We are particularly interested in the spatial structure of the signal mode a\, since one needs to know this 
spatial structure to design mode-matching in practical experiments. The integration J d 3 k can be expressed as 

/ d 3 k = J^°^^ sa ^g/2 Ji-w wnere — sin 6d9d(p represents the infinitesimal solid angle and S is the bandwidth 
of the spontaneous emission field which is in the order of the natural width T of the excited level |e). The spatial 
structure of the mode a\ is determined by the superposition coefficient /k = / {oj, fi) as a function of the solid angle 
il. Typically, the size L of the atomic ensemble (centimeter long or less) satisfies the condition TL/c <C 1, and 
oj sg L/c <C 1 {oj eg is either around GHz or zero, depending on one uses the hyperfine level or the Zeeman sublevel 
for the |s) state). In this case, Ak • = (k— fco z o) ' r i — (k/ |k| — Zq) -r^, which only depends on the direction of 
k, and becomes independent of the frequency oj. We assume that the Raman pumping field can be factorized as 
u (r, t) = u± (r) hi {z, t), where p,i {z, t) is a slowly changing function of the coordinate, and can be well approximated 
as a constant along the size of the atomic ensemble. With this condition, the superposition function / {oj, fi) is fac- 
torized as a product of the frequency part / w {oj) and the spatial part fa (17), and the spatial part /a (Tl) becomes 
independent of the emission frequency u> and the interaction time to- To have a more explicit expression of the spatial 
structure /q {fl) of the signal mode a\, we assume that the atoms are distributed by the normalized distribution 
function pdis (r) (with J p^is (r) d 3 r = 1), which is determined by the geometry of the ensemble. As we mentioned 
before, the interaction coefficient <?k normally also depends on the direction of the wave vector k in the form of a 
dipole pattern, but it varies very slowly with the angle f2 compared with other contributions in /n (f2), so we drop it 
off when considering the spatial structure of the signal mode. Under this circumstance, /n (Q) is simply expressed as 

/n(fi) =<H X (r)e- iAk - r *) {rj} (18) 

= J d 3 TU± (r) Pdis (l") e ik o z ( 1 ' cose )~ ik o sin6(x cos <p+y simp) 

To further simplify this expression, we assume a Gaussian pump beam with u± (r) = e~( x +v V r °, and a Gaussian 
form for the transverse part of the atomic distribution function with |L8| 



Pdis (r) =0, {z < -L/2, or, z > L/2) (19) 



where tq and Rq characterize the radius of the pump beam and the radius of the atomic ensemble, respectively, and 
L is the length of the ensemble. In this case, we have an analytic expression for the signal mode function 

h (O) = ^i- s e-i^' i » i, /(^) sine ffcoisin 2 t) , (20) 

where the function sine is defined as sine (x) = sin (x) jx. From this expression, we see that the signal photon 
mainly goes to the forward direction. The signal mode is inside a small cone around 8 = with A8 characterized by 

min V^o + R ll (kor Ro) , 1/ V^oL . 

We have shown that the first component of the solution (14) exactly contributes to the ideal coherent process, and 
have specified the spatial structure of the signal mode. Now we analyze the contributions of the noise described by 
the second component of Eq. (14). From Eq. (14), we see that the noise component p2Pn is expressed as a tensor 
product of the atomic density operator = J^i^i \ s )i ( s \ an d tne remaining optical density operator p v n . Besides the 
symmetric collective atomic mode s s and the signal optical mode a s , there are many other noise modes contributing 
to the atomic density operator p° and the optical density operator p^. These noise modes correlate with each other 
in a complicated way, however, these correlations do not contribute to the noise of the relevant dynamics if the two 
signal modes are not involved in the correlations, since the reduced density operator of the modes s s and a s , which 
describes the relevant dynamics, will not be influenced by these correlations after we take trace over all the noise 
modes. Nevertheless, the correlations between the good modes and the noisy mode will have influence on the relevant 
dynamics, and these contribute to the two sources of noise we have mentioned before: the spontaneous emission 
loss and the inherent mode-matching inefficiency. These two sources of noise are described quantitatively by the 
conditional probabilities p spon and p mo dc, respectively. The p sp0 n represents the probability that an atom is excited 
to the right mode s s , but the accompanying photon does not go to the signal mode a s . From the solution (14) to the 
whole density operator p/, this possibility can be expressed as 
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.Pspon 



= 1 - 



tr ((OgQpl a s s s pfslal\O a O p )) 
tr ({O a \ s s p f si\O a )") 



(21) 



where tr (• • •) represents the trace over all the remaining atomic and optical modes involved in the operator pf. 
Similarly, p mo de represents the probability that a photon is emitted to the signal mode a s , but the accompanying 
atomic excitation is not in the right mode s s , and this possibility can be expressed from the density operator pf as 



J'modc 



tr ((OqQpl a s s s pfslal\0 a 0p)) 
tr ((0 P | a s p/a!|0p}) 



(22) 



In the following, we need to calculate these two probabilities to quantify the noise magnitudes. 

From the solution (14), we can derive the probability tr ((0 o | s s pfsl\O a )) for one atomic excitation in the mode s s , 
the possibility tr ((0 P | a s pfal\O p }) for one photon in the mode a s , and the joint possibility tr ((0 a p | a s s s p/sjaj|0 a p )) 
for one atom in the mode s s and one photon in the mode a s . They are respectively given by 



tr((0 a \s s p f sl\0 a )) =p c +p 2 /N a , 
tr «0 P | a s p f al\O p )) =p c + XPc 
tr ((0 o p | a s s s p f slal\0 a 0p)) =p c + XPc/N a , 



(23) 
(24) 
(25) 



where p 2 and p c are obtained respectively from the normalization of p n and a\ ( see Eq. (16) and Eq. (17)), and the 
ratio x is defined as 



X = N a tr((O p \a s pP n al\O p })/p c , 



(26) 



with pP, the optical part of the noise component pip n as we have specified before. 

From the above equations, we see that we only need to calculate the two ratios p c /p2 and x to determine the 
noise probabilities p sp on and Pmode- As we have mentioned before, the shape u(r,t) of the pump beam is typically 
decomposed as u (r, t) = u± (r) pi (z, t) with pu (z, t) approximately independent of z along the size of the atomic 
ensemble. In this case, the two ratios p c jpi and % have much simplified expressions, which, become independent of 
the interaction details, such as the interaction time or the bandwidth of the coupling field, and depend only on some 
spatial integrations determined by the geometry of the atomic ensemble and the pump beam. For simplicity, we also 
neglect the slow variation of the coupling coefficient with the direction of k. Under these conditions, the ratios 
PcIpi and \ are given respectively by 



Pc_ 

P2 







(u ± (ri)e- iAk -") {P(} 




\ 








(u ± (ri)e- iAk 







(27) 



f dn f dir((K(r 4 )| 2 e - l ( k - k ')- r A 

« (r 4 )e lAk - r ') {ri} («± (r<) 



(28) 



{r 4 } 



1. 



Similar to the case for calculating the signal mode structure /q (f2), we also assume a Gaussian pump beam and a 
Gaussian form for the transverse atomic distribution function as is shown by Eq. (19). In this case, we have the 
following analytic expressions for these two ratios 
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Pc ^ rj (rl + 2Rl) 
P2 2{rl+Rlf 



x sine (k L sin 2 - ) J- , (29) 



X ■ 



sin( 



OdO \ e -hklrlRl^6/(rl+Rl) ^ / ^ ^ 



1 -2 



^ (rg + 2i?2) 

X £ £ Sin Sin O'dOdO' j i- jf ^ e^o^ (-" 2 »'-2 -n 9 sin 9' cos V )/ (r 2 +2i? 2 ) 



(30) 



xe 



■±/c 2 r- 2 i? 2 (sin 2 9+sin 2 9')(r 2 +_R 2 ) 



sine 



i t I ■ 2 u ■ 2 ° 

k L sin sin - 

2 2 



sine [ k a L sin 2 ^ I sine [ fc ^ sin 2 tt 



1. 



In writing Eq. (29), we have assumed koL >• 1 and have neglected the term which is about l/k^L times smaller. 

From the two ratios p c /p2 and \, the spontaneous emission loss and the inherent mode matching inefficiency are 
directly written as 



Pspon ^ (1 + N a p c /p 2 ) 1 



(31) 



Pmodo - xl (1 + x) • (32) 

The approximations in Eqs. (31) and (32) are valid under the condition N a 3> 1 (we neglected the terms which are 
1/N a times smaller). To minimize the two sources of noise p spon and p mo dc, it is better to have a large N a p c /p2 and 
a small %. We will discuss the details in the following section on how to minimize these two noise by using different 
kinds of mode matching methods. 

At the end of this section, we would like to make a brief comparison with the single atom case. For the case of a 
single atom, if the position of this atom is also fluctuating in the space, the above calculation is still valid but with 
N a = 1. From Eqs. (21)-(26), we see that in the case of 7V a = 1, we always have p mo do = 0. Thus, for the case of a 
single atom, the spontaneous emission loss is the only source of noise, as has been shown intuitively in Fig. 4. However, 
in this case, the spontaneous emission inefficiency p spon , which can be approximated by p spon ~ 1 — (1 + x) Pc/P2 since 
Pc ^ P2i becomes much larger. If we compare the efficiency (defined as 1 — p spon ) from the spontaneous emission 
noise between the single-atom case and the atomic-ensemble case, we see that this efficiency is increased by about a 
factor of N a for the atomic ensemble. This is what we called the collective enhancement of the signal-to-noise ratio 
for the many-atom ensemble. 



V. MODE MATCHING METHODS 



The first mode-matching method is to choose the signal mode as the mode for detection, with the mode function 
exactly given by the inherent spatial structure (20) of the signal light. The mode function (20) is similar to a Gaussian 
function, especially when k a r^ > L and fc i?ci > L. (It would be an exact Gaussian function if the axial distribution 
of pais (r) is Gaussian) . One can couple a Gaussian mode into a single- mode optical fibre with a good technical 
mode-matching efficiency. For this case with an exact mode matching, the two noise probabilities p sp0 n and p mo dc are 
exactly given by Eqs. (29)-(32). We now go to the calculations of these two noise probabilities. For this purpose, we 
need to understand first how to get a large N a p c /p2 and a small \ to minimize the noise p S pon and p m odo- 

First, we give an estimation of the ratio N a p c /p2- It has a definite physical meaning. From Eq. (29), we see that 
the integration function is significantly different to zero only when the integration variable 9 falls into a small cone 
with < 6 < Of, where 9f is estimated by Of ~ min (l/^/fcoL, l/fcor e //) with r e // = r -Ro/\Ao + ^o- Thus, the 
result of the integration can be estimated by 2 j. We also know that the total atom number can be approximated 
by N a ~ n a irRlL, where n a denotes the average atom number density. From this, we see that the ratio N a p c /p2 is 
estimated by 

NaPc . (I . 2 1 \2 r \ 

— — ~ mm ^-n a X r ef f, ^n a \L j . (33) 
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If we assume that the two bounds for N a p c /p2 are comparable, which means that l/^/k^L ~ l/kor e ff, or the Fresnel 
number defined as Fr = nrl^/XoL is in the order of 1, we have N a p c /p 2 ~ n a \\Lj (47r) ~ d /(47r), where the 
on-resonance optical depth d a is simply defined as d = n a \\L. Therefore, to have a small spontaneous emission 
inefficiency p S pon, we need to use an optically dense ensemble with d Q 3> 1 (the off- resonance optical depth can still 
be much smaller than 1 due to the large detuning). 

Next, we consider how to minimize the ratio x by choosing appropriate interaction configurations. From Eq. (30), 
we see that x increase with the ratio Rq/tq. This can be intuitively understood as follows: the noise characterized 
by Pmode and x comes from the density fluctuation of the atomic ensemble. With a larger Rq, the atoms have more 
space to move around, and one thus has a relatively larger density function and a larger noise ratio x ]l9| - Thus, to 
minimize x, we choose a configuration with ro S> Ro in the following, which means, the pumping laser is shined on 
all the atoms with a broad cross section. In this case, the effective radius r e ff ~ Rq. 

Now we would like to calculate the noise probabilities p sp0 n and Pmode m a more accurate way in the case of exact 
mode matching. This requires us to carry out the complicate integrations in Eqs. (29) and (30), which is only possible 
by using numerical methods. The spontaneous emission inefficiency p spon depends on the optical depth d — n a \\L. 
We have numerically calculated the value of p S pon versus the optical depth d under different geometries of the atomic 
ensemble. The geometry of the ensemble is described by the Fresnel number Fr = ttRq/XqL. The calculation result 
is shown in Fig. 5. From the figure, we see that the spontaneous emission inefficiency is insensitive to the geometry 
of the ensemble. The two curves with Fr = 1 and Fr = 10 basically overlaps. With a much smaller Fr — 0.1, the 
inefficiency p S pon increases by about a factor of 2. On the other hand, p S pon is sensitive to the optical depth d . It 
decrease with l/d approximately linearly, which confirms the rough estimation from the above. From the curves, we 
can write p S p 0n approximately asp spon = l/(l + d /26) for the Fresnel number F r from 1 to 10. In the above numerical 
calculation, we assumed some typical values for the parameters with Ao ~ 0.8/im and L ~ 1cm, but actually the result 
is very insensitive to the parameter kgL as long as it is still much larger than 1. For an atomic cell with a typical 
number density n a ~ 10 12 /cm 3 and length L ~ 1cm, the optical depth d Q ~ 6.4 x 10 3 , and we have a spontaneous 
emission inefficiency p sp0 n ~ 0.4%, which is very small. 

0.35 1 1 1 1 1 1 1 1 

0.3 \ 

0.25 -1 

0.2-1 
Pspon . ' 
0.15 1 



0.05 



2000 4000 6000 8000 10000 12000 14000 

d 

FIG. 5. The spontaneous emission inefficiency p S pon versus the on-resonance optical depth d of the ensemble under different 
geometries characterized by the Fresnel number Fr. We have Fr — 10 for the solid curve, Fr = 1 for the dotted curve, and 
Fr = 0.1 for the dashed curve. 

The inherent mode matching inefficiency p mo dc is determined by the geometry of the ensemble, and docs not depend 
on the value of the optical depth. We have numerically calculated the value of p mo de- It turns out that p mo dc is actually 
also insensitive to the geometric parameters koL and Fr as long as they are still in the typical region. For instance, 
we have p m odc — 24% for Fr = 1, p mo dc — 25% for Fr — 10, and p mo de — 23% for Fr = 0.1. In each case, we change 
the length L from 0.3cm to 3cm (with Ao ~ 0.8pm), and the value of p mo de basically does not change at all. Note 
that the typical inherent mode matching inefficiency p m ode — 24% is quite large, and we can not efficiently reduce 
it by controlling the geometry of the ensemble. As we have mentioned before, the mode mismatching noise has a 
more significant influence on the efficiency of the application schemes than the spontaneous emission noise. The large 
inherent mode matching inefficiency is a disadvantage of the exact mode-matching method. 

To reduce the inefficiency p mo dc, we consider the following simple improvement to the exact mode matching method. 
After the light-atomic-ensemble interaction, the scattered light is focused by a lens, and we apply an aperture on the 
focal plane of the lens to collect the signal light only from a small cone in the forward direction with < 6 < 9d , where 
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the detection angle 9d can be controlled by the size of the aperture. We call this kind of mode matching the filtered 
exact mode matching. By the exact mode matching, we detect the photon scattered to the signal mode; and by the 
filtered exact mode matching, we only detect the scattered photon which is in the signal mode and at the same time 
lies in the filtering cone with < 9 < 9d- Due to the additional restriction, we now have a less success probability to 
register the photon. The spontaneous emission inefficiency Pspon 

should increase, but through this scarification, we 
can significantly reduce the noise p mo de, which means, if we register the photon with this additional restriction, we 
are more confident that the accompanying atomic excitation will go to the collective atomic mode s s . To calculate 
p S p 0n and p m ode for the filtered exact mode matching, we can still use Eqs. (21)-(26), but the signal mode a s should 
be confined inside the small cone with < 9 < 9d- Inside this filtering cone, the signal mode has the same mode 
structure as given by Eq. (20). With this modified signal mode, we find that p mo de is still given by Eqs. (32) and 
(30), but all the integrations of 9 and 9' in Eq. (30) can only be taken from to the filtering angle 9d instead of from 
to 7r. The spontaneous emission inefficiency p sp0 n now has the following expression 

Papon ^ 1 - -T- T : , (34) 

N a Pc + P2 

where we have assumed N a 3> 1 for the approximation, and the symbol p c \e<e D means that the integration of 9 in p c 
can only be taken from to Od- From Eq. (34), we need to calculate two ratios p c /p2 and p c \8<e D /p2 for p sp on- The 
ratio p c /p2 is exactly given by Eq. (29), and the ratio p c \e<e D /P2 has the same form as Eq. (29), but the integration 
of 9 is only from to Od- 

We have numerically calculate the noise probabilities p sp on and p mo de versus the filtering angle 9d, and the result 
is shown in Fig. 6 with the Fresnel number Fr — 1 and the optical depth d„ ~ 1.9 x 10 3 . From the figure, we see 
that by decreasing the filtering angle 9d, we can significantly reduce the noise p mo de- The cost is that at the same 
time the inefficiency p sp on significantly increases. However, as we have mentioned before, though both of the noise 
p S pon and Pmode are correctable in the application schemes in Refs. |lj,[l3|], the mode mismatching noise p mo de has a 
more severe influence on the final efficiency of the schemes than the spontaneous emission noise p sp on- Thus, for these 
applications, it is worthy to choose an appropriate filtering angle 9d with p mo dc significantly smaller than p sp on- For 
instance, if we choose the angle Od — 0.002, the mode mismatching noise p mo dc — 0.9%, which is basically negligible 
compared with other sources of noise. At the same time, p sp0 n — 66.6%, which seems to be quite large. However, to 



overcome this noise, we need only to increase the repetitions of the entanglement generation scheme in Refs. [|12 13 
by another factor of 3 (given by 1/ (1 — p sp on)), which is actually only a moderate cost. This choice would be much 
better than the case without any filtering, where we have an inherent mode mismatching noise p mo de — 24%. For the 
purpose of guiding the choice of the best experimental configuration, we list several important values of p sp0 n and 
Pmode versus the angle 9d in the following table 
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FIG. 6. The spontaneous emission inefficiency p spon (solid curve) and the inherent mode mismatching probability p mo d c 
(dashed curve) versus the filtering angle 6d with the Fresnel number Fr — 1. 

We have also calculated the noise probabilities p sp on and p mo dc versus the angle Oo under different geometries of 
the ensemble characterized by the Fresnel number Fr. For Fr = 10 and Fr = 0.1, the results are shown in Fig. 7. In 
all the calculations, we assumed the same optical depth d„~1.9x 10 3 and the cell length L ~ 1cm. The qualitative 
properties of the curves for different Fresnel numbers are quite similar, but one needs to shift the appropriate angle 
9r>- With a large Fresnel number Fr, one needs to choose a smaller angle 9n- The appropriate filtering angle 9d 
changes by a factor of 2-3 if the Fresnel number changes by a factor of 10. 
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FIG. 7. The spontaneous emission inefficiency p sp0 n and the inherent mode mismatching probability p mo de versus the filtering 
angle 9d with the Fresnel number Fr = 10 (the two dashed curves) and Fr = 0.1 (the two solid curves). The spontaneous 
emission inefficiency increases with the decrease of the angle 6d while the mode mismatching noise decreases. 

Finally, we describe a mode matching method which is simplest for implementation. For this method, we still use 
an aperture to select the scattered light from the small cone with < 9 < 9q. But now we do not use any technical 
mode matching to choose only the signal mode with the right mode structure /n (Q) for detection. Instead, we detect 
the photon in all the modes which lie in the small cone with < 9 < 9o. We call this the simple filtering method. This 
seems to be a bad mode-matching method, since one does not make any distinction between the modes in this small 
cone. If the registered photon does not come from the right mode a s , the accompanying atomic excitation will be 
most probably not in the collective mode s s , and one thus has a significantly larger inherent mode mismatching noise 
Pmade- However, the observation here is that if 9jj is sufficiently small, basically only the signal mode exists in this 
small cone, and all the other modes make negligible contributions. The important question is how small the angle 9d 
should be to guarantee a small p mo dc- Note that for the simple filtering method, the spontaneous emission inefficiency 
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Pspon can be calculated in the same way as the filtered exact mode matching method by the use of Eq. (34). But we 
have a different expression for p mo do, which in this case is defined as the probability that the accompanying atomic 
excitation does not go to the collective mode s s when one photon (from any optical modes) is detected in the small 
cone with < 9 < 6d ■ This possibility can be expressed as 

, Pc\e<0 D + (1/jVo) tr (p 2 pn) \e<e D , Q rX 

Pc\e<0 D + tr (P2Pn) \0<0 D 

where the trace tr (• • •) is over all the atomic and the optical modes, and the symbol \g<g D has the same meaning that 
the integrations of 9 in p c and tr (j>2Pn) are only from to 9d- In the limit of -ZV a ^> 1, p mo de can still be written into 
the form p modc ~ x«/ (! + Xs)> with 



Pc\6<e D 



Xs + 1 tr (p2Pn) \0<0d +Pc\0<9 D 



1 



(1 — cos9d) 



sm.6d6 







sin 2 



x sine ^k L sin 2 | . (36) 

In writing Eq. (36), we have assumed the same kind of atomic distribution function pdis (r) and the pump mode 
function u± (r) as we have specified before. We can use Eqs. (36) to numerically calculate the mode mismatching 
probability p mo do versus the filtering angle 9]j for the simple filtering method. The result, together with the curve 
for p S p 0n versus 9d, is shown in Fig. 8 for the Fresnel number Fr = 1. In the calculation, we assumed the same 
optical depth d ~ 1.9 x 10 3 . From the figure, we see that if 9d is large, the mode mismatching probability p mo do is 
very large, but as 9b decreases, p mo do can still tend to zero. This confirms our intuitive observation. With the same 
parameters, for the simple filtering method, we should further decrease the filtering angle 9u by a factor of 2-3 to get 
the optimal configuration compared with the filtered exact mode matching method. Of course, due to the decrease 
of the optimal 9]j, the corresponding spontaneous emission inefficiency p sp0 n significantly increases. To guide the 
experimental choice of the optimal 9jj , we also list some important values of j3 spon and p mo do f° r the simple filtering 
method 
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FIG. 8. The spontaneous emission inefficiency p spon (dashed curve) and the inherent mode mismatching probability p mo d c 
(solid curve) versus the filtering angle 6d with the Fresnel number Fr — 1 for the simple filtering method. 



For the simple filtering method, we also calculated the noise probabilities p S pon and p mo dc under different geometries 
of the ensemble. For the Fresnel number Fr — 10 and Fr = 0.1, the results are shown in Fig. 9. We assumed the 
same optical depth and the cell length in these calculations. Since the qualitative picture revealed by this figure is 
quite similar to the case of the filtered exact mode matching, we do not need to discuss its properties in details. 
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FIG. 9. The spontaneous emission inefficiency p ap on and the inherent mode mismatching probability p mo do versus the filtering 
angle 6d with the Fresnel number Fr = 10 (the two dashed curves) and Fr = 0.1 (the two solid curves) for the simple filtering 
method. 



For experimental realizations, the simple filtering method can be much easier than other mode matching methods. 
The calculation here shows that the price we need to pay for this simplification is to choose a smaller filtering angle with 
a significantly increased spontaneous emission inefficiency. Note that for the application schemes in Refs. jl2],[l3), this 
price seems to be still acceptable. For instance, if we choose 8rj — 0.001, the mode mismatching noise p mo dc — 4.8% 



has already been small. In this case, the additional repetitions of the entanglement generation scheme in Refs. [|12 13 
is given by the factor 1/ (1 — p sp0 n) — 10, which is not too much. 



VI. SUMMARY 



In summary, we have developed a theory in the weak pumping limit to describe the three-dimensional effects in the 
interaction between the free-space light and the many-atom ensemble. The calculations demonstrate some interesting 
results which are not know from the other approaches: firstly, it shows that the signal light has an inherent spatial 
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mode structure, which is determined together by the geometry of the ensemble and the mode structure of the pump 
beam. Secondly, it reveals that there will be two sources of noise during the light-atom interaction. One is the 
spontaneous emission inefficiency, which is inversely proportional to the on-resonance optical depth of the ensemble; 
and the other is the inherent mode mismatching noise, which arises from the density fluctuation of the ensemble, and 
should be fully determined by the geometry of the ensemble for room-temperature atomic cells. It turns out that 
the inherent mode mismatching noise is quite large if one collects the signal light through the exact mode matching 
method, and one cannot efficiently reduce this noise by optimizing the geometry of the ensemble since in the typical 
parameter region, this noise is insensitive to the cell geometry. Finally, we show an effective way to reduce the inherent 
mode mismatching noise by adding an aperture to select the light only from a small emission cone. By this mean, 
we can efficiently reduce the mode mismatching noise at the price of increasing the spontaneous emission inefficiency. 
It is worthy to do this since the spontaneous emission inefficiency is far less important for some application schemes. 
There are two methods for this purpose: the filtered exact mode matching method and the simple filtering method. 
Both methods can reduce the inherent mode matching noise if one chooses an appropriate filtering angle. 

The calculations in this paper show that there is a large collective enhancement of the signal-to-noise ratio if 
one compares the many-atom ensemble with the single-atom case, which confirms the predication from some simple 
theoretical models, and demonstrates the validity of this important observation for the complicated realistic situations. 

At the end of this paper, we should mention that the calculations here also show that for any mode matching 
method, we can not make the spontaneous emission inefficiency Pspon 

and the inherent mode mismatching noise p mo de 
both negligible for the room-temperature atomic cell even if the cell has a large on-resonance optical depth. This 
result does not have much influence on the application schemes in Refs. |12 since they are inherently robust to 
the noise p sp0 n and p mo dc- As a consequence of this inherent robustness, a non-negligible but not-huge noise only 
has a moderate influence on the efficiency of the schemes. However, for any of the application schemes of atomic 
ensembles which are not inherently robust, such as the quantum light memory scheme [0J§[?J^] or the continuous 
variable teleportation scheme a non- negligible noise will decrease the scheme fidelity. The calculations here do 

not directly apply to these schemes, since they work out of the perturbation region. However, it indeed raises the 
interesting question whether there is some theoretical limit on the best achievable fidelity of these schemes if one uses 
free-space atomic cells for their implementation. 
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cell and for the pump beam, but they are insensitive to the detailed form of the atomic distribution function. Here, for 
simplicity of the calculation, we assume a Gaussian form for the atomic transverse distribution function. But we expect 
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[19] Note that this is a valid picture for room-temperature atomic cells. However, for ultracold atoms, each atom will be 
localized in a small space, and the density fluctuation basically will not depend on the size and geometry of the ensemble. 
Thus, for cold atoms, one may have a different interaction picture and a smaller noise ratio. 
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